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We present a preliminary analysis of the heavy-heavy spectrum and heavy-light decay constants in full QCD, 
using a tadpole-improved SW quark action and an RG-improved gauge action on a 16 3 x 32 lattice with four sea 
quark masses corresponding to m n /m p w 0.8, 0.75, 0.7, 0.6 and a -1 ~ 1.3 GeV. We focus particularly on the effect 
of sea quarks on these observables. 



1. Introduction 

Despite the difficulties that heavy quarks 
present for study on the lattice, in two areas at 
least, they provide an excellent laboratory for ex- 
amining effects of sea quarks, namely the hyper- 
fine splitting of hcavy-onia and heavy-light decay 
constants. 

The calculation of the hyperfine splitting has 
been extensively pursued in quenched QCD0. 
In the charm sector, one can compare results 
with the experimental value Mj/^ — M nc = 
118(2)MeV. Relativistic actions underestimate 
this splitting, giving a result of 70-80 MeV. Non- 
relativistic actions yield similar values. This dif- 
ference has been argued to arise from quench- 
ing, due to a smaller value of the strong coupling 
constant at short distances in quenched QCD[Q. 
Since the strong coupling constant, as estimated 
from the Coulomb term of the static potential, 
has been seen to increase by about 10% by the 
introduction of sea quarks one may hope to 
observe an effect also in the hyperfine splitting. 

A similar reasoning suggests that heavy-light 
decay constants are underestimated without sea 
quarks. Recent results from the MILC H and 
NRQCD H Collaborations do show that the de- 
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Table 1 

Run parameters for this simulation with eventual 

final statistics in parentheses. 

K aea mps/my a- 1 [GeV] N cfg (HH,ffl7T 

0.1375 0.8048(9) 0.9653(65) 179, 322 (700) 

0.1390 0.751(1) 1.0205(80) 215, 311 (700) 

0.1400 0.688(1) 1.0889(72) 270, 229 (700) 

0.1410 0.586(3) 1.1612(87) 169, 224 (500) 

cay constants increase with decreasing sea quark 
mass, although an accurate measurement still has 
to be carried out. 

In this report we present preliminary heavy 
quark results calculated on two-flavor full QCD 
configurations being generated by the CP-PACS 
Collaboration^, with emphasis on search for sea 
quark effects. 

2. Computational Details 

Our analysis is carried out on a 16 3 x 32 lat- 
tice for four sea quark masses corresponding to 
mps/mv ~ 0.8-0.6. Configurations are gen- 
erated with an RG-improved gauge action at 
j3 = 1.95 and the SW quark action with the 
clover coefficient c sw = P~ 3 / 4 = 1.53 where 
P = 1 — 0.8412//3 is the one-loop value of plaque- 
tte. The lattice spacing, set either by string ten- 
sion or m p , takes the value a -1 w 1.3 GeV. Pa- 
rameters relevant for the present study are listed 
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Figure 1. Heavy hyperfine splitting using kinetic 
mass definition. Scale is set by string tension 
measured for each sea quark mass. 

in Table @. 

Hadron correlators are computed using local 
and exponentially smeared sources in Coulomb 
gauge. For study of the heavy- heavy spectrum we 
choose values of K corresponding to araq — 0.68, 
0.75, 0.83 and 0.9 in the naive pole mass definition 
arriQ = (1/ K — l/K c )/2. In conjunction with the 
above estimate for a -1 , our present study there- 
fore explores the region of charm quark. 

For the heavy-light decay constants a set of 
lighter heavy quark masses are chosen, deter- 
mined from a plot of the kinetic mass of heavy- 
light meson against 1/K. The light quark has 
a mass tuned to approximately the strange quark 
mass as estimated from m~f /m^ — 0.688 accord- 
ing to the Gell-Mann-Okubo formula. 

3. Results 

3.1. Hyperfine splitting 

We calculate the hyperfine splitting from the 
difference of pole masses. With the use of the 
SW action lattice discretization effects associated 
with large quark masses are expected to be small 
for this quantity. The pole mass itself, how- 
ever, does suffer from large discretization effects, 
causing problems in calculating the ground state 
mass. Several possible schemes for measuring the 
ground state mass have been proposed. We corn- 



Figure 2. Heavy hyperfine splitting using HQET 
mass definition. 

pare two : 

(i) the kinetic mass, derived from the dispersion 
relation, 

(ii) an "HQET" mass, suggested by Bernard et 
aZ.[p|, defined as 

a m HQET = am pole + #Q (<»"*«„ ~ am pole) '. W 

with Nq the number of heavy quarks. We 
use tree-level quark mass definitions tadpole- 
improved by a factor uo = P 1 / 4 . 

We show our results in Figs. [I] and | The 
scale is set by the string tension determined 
for each sea quark mass to absorb change of 
scale. As can be seen, the two definitions draw 
different conclusions for the hyperfine splitting. 
The kinetic mass gives results which indicate an 
almost negligible change between the heaviest 
sea quark mass (~ 3m strange) and the lightest 
(~ m strange/ '2) . The values are roughly consis- 
tent with those obtained from quenched stud- 
ies. Using the HQET mass definition there is 
a significant difference between the heaviest and 
lightest sea quark masses, but the pseudoscalar 
masses are now much smaller, indicating that at 
the charm quark mass, the hyperfine splitting is 
smaller, approximately 55 MeV. 

3.2. Heavy-light decay constants 

In Fig. H we plot fpsV&Tps, using the KLM 
normalization for quark fields. The axial vector 
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Figure 3. fps-\/rnps versus l/mps with a valence 
strange quark. For clarity only the heaviest and 
lightest sea quark masses and their central fits are 
shown. 

renormalization fatcor is taken from Ref. H cal- 
culated for massless quark at one-loop perturba- 
tion theory. As the heavy quark masses taken are 
smaller than those for the hyperfine splitting mea- 
surements, the difference between the two mass 
definitions is negligible. We adopt the HQET 
mass definition which has smaller statistical error. 
The scale is set by the string tension as before. 

The results for the heaviest and lightest sea 
quark masses are quite similar in magnitude. 
There is a slight upward change in the slope, how- 
ever, for lighter sea quark. An extrapolation to 
obtain fn a , for which we adopt an expansion lin- 
ear in l/mps, yields a shift between the heaviest 
and lightest sea quark masses of 2-3 %. 

4. Conclusions 

At present the ambiguity that exists in the 
ground state meson mass precludes any quanti- 
tative statement about sea quark effects in the 
hyperfine splitting. We emphasize that this is 
not symptomatic of the dynamical configurations 
employed but is a general feature of coarse lattice 
spacings. This problem will be alleviated in the 
set of configurations with a finer lattice spacing 
a" 1 ~ 2.5 GeV generated by CP-PACS, which we 
plan to analyze next. Another possible avenue of 
approach would be to reduce discretization effects 



through introduction of further terms in the SW 
action as suggested in Refs. |]lo| , p| . 

The heavy-light decay constants around the 
charm mass region examined here show only 
a slight shift upwards as the sea quark mass 
is reduced. Quenched chiral perturbation the- 
ory [Il2 13 1 suggests, however, that the shift in- 



creases with the heavy quark mass, and may reach 
a sizable level in the region of the b quark. Stud- 
ies with very heavy (or static) quarks may then 
provide the best region for understanding the sys- 
tematic error due to quenching. 
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